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ABSTRACT 

We present high-precision radial velocity observations of WASP- 17 throughout the transit of its 
close-in giant planet, using the MIKE spectrograph on the 6.5m Magellan Telescope at Las Campanas 
Observatory. By modeling the Rossiter-McLaughlin effect, we find the sky-projected spin-orbit angle 
to be A = 167.4 ± 11.2 deg. This independently confirms the previous finding that WASP-17b is on a 
retrograde orbit, suggesting it underwent migration via a mechanism other than just the gravitational 
interaction between the planet and the disk. Interestingly, our result for A differs by 45 ± 13 deg from 
the previously announced value, and we also find that the spectroscopic transit occurs 15 ± 5 min 
earlier than expected, based on the published ephemeris. The discrepancy in the ephemeris highlights 
the need for contemporaneous spectroscopic and photometric transit observations whenever possible. 
Subject headings: planetary systems: formation - stars: individual (WASP-17) - techniques: radial 
velocities 



1. INTRODUCTION 

Knowledge of the angle ip between a star's spin vector 
and a planet's orbital angular momentum vector con- 
tains information about the planetary formation process. 
This is especially so for Hot Jupiters, which are thought 
to have "migrat ed" inward from large orbital distances 
(|Lin et al.l [l996) . One possible migration mechanism in- 
volves the gravitational interaction between t he planet 
and t he gas in the pr otoplanetary disk (e.g. iLin et all 
119961: llda fc Lin! 12004ft . Such interactions would pre- 
serve a prograde orbit, and could not account for Hot 
Jupiters with retrograde orbits (ip > 90°) unless the pro- 
toplanetary disk was initially misaligned with the star 
(jBate et al.1 I2010I ). Other migration mechanisms, in- 
volving planet- planet scattering and Kozai cycles, (e.g. 
iNagasawa eFaXl [2008) can produce retrograde orbits. 
Therefore by measuring tjj one gains insight into the mi- 
gration history of a particular planetar y system. 

The Rossiter-McLaughlin effect (jRossiterl Il924t 
IMcLaughlinl fl92l occurs when part of the rotating 
stellar photosphere is eclipsed by a companion star 
or planet. This removes a velocity component of the 
observed rotationally- broadened line profiles, causing 
a pattern of anomalous radial velocities to be observed 
throughout the eclipse. Although this effect does not 
reveal it does give a measure of A, the angle between 
the sky projection of the stellar spin vector and the 
planetary orbital angular momentum. The Rossiter- 
McLaughlin effect has now been used to measure A for 
28 transiting plan et systems (see lTriaud et al.l ((2010h or 
I Winn et allMlOh L 

The tran siting Hot Jupiter WA SP- 17b, which was dis- 
covered by lAnderson et all (|2010h . is reported to have a 
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planetary radius of Rp — 1.74 Rj and a host star ro- 
tational velocity of wsini = 9km s . The deep transit 
and fast stellar rotation make it a propitious target for 
R ossiter-McLaugh l in effe ct observations. 

lAnderson et al.l (|2010T l reported that the planet ap- 
peared to be in a retrograde orbit, based on three radial- 
velocity measurements that were obtained during tran- 
sits. Subsequently, the same group announced further 
radial velocity measurements of W ASP-17, with denser 
time sampling (jTriaud et al.l I2010T ). Those data gave 
much stronger support to the claim of a retrograde or- 
bit, and indeed g ave a value for A nearly identical to the 
value derived by lAnderson et al.l (|2010f) . In this Letter 
we present new observations and an independent confir- 
mation of the likely retrograde orbit of the WASP-17 sys- 
tem, though with a significantly different derived value 
for A. 

2. OBSERVATIONS 

On the night of 2010 May 11, we obtained high- 
resolution spectra of WASP-17 using the Magellan In- 
amori Kyocera Echelle (MIKE) spectrograph on the 6.5m 
Magellan II (Clay) telescope at Las Campanas Obser- 
vatory. MIKE was used in its standard configuration 
with a 0.35" slit and fast readout mode, which deliv- 
ered 32 red echelle orders with wavelengths from 5000- 
9500 A and with a resolving power of R ~ 48, 000. Spec- 
tra were taken continuously for 9 hours, with an interrup- 
tion of approximately 30 minutes when the star crossed 
the meridian near the zenith. A total of 33 spectra of 
WASP-17 were obtained, with 600 s exposures. ThAr 
calibration frames bracketed each spectrum. Observing 
conditions were good, although some high, thin clouds 
were present in the latter half of the night. The moon 
was down throughout the WASP-17 observations, and 
the seeing was 0.5-0.7". 

The spectra were reduced using the Carnegie MIKE 
pipeline developed by Dan Kelson. Radial velocities were 
determined via cross correlation with respect to a tem- 
plate spectrum. We tested different choices for the tem- 
plate spectrum, including a radial velocity standard ob- 
served on the same night and various individual WASP- 
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17 spectra. The lowest inter-order scatter was obtained 
when using the highest signal-to- noise WASP- 17 spec- 
trum as the template. This spectrum was also obtained 
at low airmass (1.005). For the cross-correlation analysis 
we selected 13 echelle orders that were free of obvious 
telluric lines and had a signal-to-noise ratio exceeding 
100 pixel^ 1 . Cross-correlation was performed using the 
IRAJ0 task fxcor. We used the O2 absorption features 
from 6870-6900 A to provide a constant reference for our 
wavelength solutions for each spectrum. These absorp- 
tion bands have been show n to be stable at the 5 m s _1 
level over short timescales ()Figueira et al.|[20Toh . 

MIKE lacks an atmospheric dispersion corrector, and 
this resulted in an airmass-dependent systematic trend 
in the radial velocity measurements, with an amplitude 
that increased for bluer echelle orders. The continuous 
nine-hour monitoring allowed us to track this systematic 
trend and account for it during the model-fitting pro- 
cedure, as described in Section [3J The order-averaged, 
airmass-corrected radial velocities are given in Table 1, 
with uncertainties taken to be the standard deviation of 
the mean of the results of all 13 orders. 

3. ANALYSIS 

Our model for the radial velocity data has the form: 

VcalcnC*) = Krb(*) + ^RM(*) + C0+Cin+(C2 + C 3 n)X, (1) 

where V ca i Ctn (t) is the calculated radial velocity at time 
t in echelle order n (ranging from 1 to 13), V OT b is the 
radial velocity due to the star's orbital motion (assumed 
to be circular), Vrm is the transit-specific "anomalous 
velocity" due to the Rossiter-McLaughlin effect, and 
{co, 01,02,03} are constants specifying the offset between 
the relative barycentric velocity of WASP- 17 and the ar- 
bitrary template spectrum that was used to calculate ra- 
dial velocities. To account for the wavelength-dependent 
airmass trend mentioned in the previous section, the off- 
set was allowed to be a linear function of both order 
number n (effectively a wavelength index) and airmass 
X. 

Many approaches have been taken to model the 
Rossiter-McLaug hlin effect, such as th e pixellated photo- 
spheric mode l oflQueloz et al.l jj2000D , the first-momen t 
approach of lOhta et al.l (|2005f ) and iGimenezj (120061) . 
the forward-modeling appro ach oflWinn et all (120051) or 
the spectral line fitting of JAlbrecht et al.l <|2009f ) and 
iCollier Cameron et all (|2010» . In our case the most ap- 
plicable and c onven ient model is the analytic formula of 
iHirano et all (|20Tol ). 



V RM {t) = -5(t) v p (t) 
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where S(t) is the loss of light during the transit, v p (t) is 
the mean radial velocity of the small portion of the pho- 
tosphere that is hidden by the planet, vsini is the stellar 
projected rotation rate, and vq is the velocity width of the 

4 IRAF is distributed by the National Optical Astronomy Ob- 
servatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 
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Fig. 1. — The spectroscopic transit of WASP-17b. Top. — 
Radial velocity variation observed on the night of 2010 May 11- 
12, along with the best-fitting model. The small gray points are 
the order-by-order radial velocities that were fitted. The larger 
blue points are averages of the results from all orders. Bottom. — 
Residuals between the data and the best-fitting model. 

spectral lines that would be seen from a small portion of 
the photosphere (i.e. due to mechanisms other than rota- 
tion) . This formula relates the phase of the transit to the 
radial velocity anomaly derived from cross-correlation. 

To calculate v p (t) we assumed that the stellar photo- 
sphere rotates uniformly with an angle A between the 
sky projections of the spin vector and the o r bital an- 
gular momentum vect or (see, e.g.. lOhta et al.l (|2005l ) or 
IGaudi fc Winnl (120071) 1 We took vsini to be a free pa- 
rameter, and vq was constrained by the condition that 
the quadrature sum of vq and vsini s hould equal the 
obser ved total linewidth of 11 km s _1 (|Anderson et all 
120101) . 

To calculate S(t), we assumed a linear limb-darkening 
law with a fixed coefficient of 0.7, and used the analytic 
formulas o f iMandel fc Ago! (|2002D . as implemented by 
iPall J2008). The parameters of the photometric model 
were the planet-to-star radius ratio R p /R- kl orbital in- 
clination i, and normalized stellar radius R*/a (where a 
is the orbital distance). Some of the transit character- 
istics are more tightly constrained by previous observa- 
tions of photometric transits than by the radial velocity 
data presented here. Hence, we used Gaussian priors for 
the transit depth (i? p /i?*) 2 , total duration t\v — ti, par- 
tial duration fry - tin, orbital period P, and time of tran- 
sit T c based on the results presented by lAnderson et all 
(2010) (see their Table 4, case 3). We also used those 
results to set a Gaussian prior for K, the orbital velocity 
semi-amplitude, since our coverage of the spectroscopic 
orbit is sparse. 

All together there are 11 adjustable parameters in our 
model fit, of which 7 are mainly determined by our new 
radial velocity data: A, vsini, cq, c±, 02, 03, and the 
time of conjunction T c . The other 4 parameters, R p /R±, 
i, R- k /a 1 and K, are controlled mainly by the priors. 
We determined the best-fitting parameter values and 
their 68.3% confidence limits with a Monte Carlo Markov 
Chain (MCMC) algorithm that w e have described else- 
where (see, e.g.. iWirm et~all(|2007f )). The likelihood was 
taken to be exp(— x 2 /2) with 
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where v b s {ti,n) is the radial velocity measured at time 
ti in echelle order n, and v ca \ c (ti,n) are the calculated 
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Fig. 2. — Illustration of airmass correction. The radial ve- 
locities for cchclle orders 61, 59, and 57 (from top to bottom), along 
with the best-fitting model that includes both an airmass correc- 
tion and the Rossiter-McLaughlin effect. Central wavelengths of 
the orders are indicated on the plot. The lower panel shows the 
time dependence of the airmass. 

radial velocities. We found that the choice a — 92 m s -1 
gave x 2 — Nd f for the best-fitting model, and is also 
approximately equal to the scatter observed in the resid- 
uals to the best-fitting model, so we used this value in 
the MCMC analysis. 

Table [5] gives the resulting best-fit parameter values. 
Figure [T] shows the airmass-corrected radial velocities as 
a function of time, along with the best-fitting model and 
the residuals. Figure [5] illustrates the airmass correction 
for three separate echelle orders. Figure [3] displays the 
marginalized, joint a posteriori probability distributions 
involving A and the other model parameters. 

The radial velocity data clearly exhibit an "upside- 
down" Rossiter-McLaughlin effect: an anomalous 
blueshift for the first half of the transit, followed by 
an anomalous redshift for the second half. Our re- 
sult for the projected stellar rotation velocity, v sin = 
8.61 ± 0.45 km s " 1 , is in a greement with the results of 
lAnderson et~aT1 (|2010fl and ITriaud et al.l (|2010D . How- 
ever, our finding of A = 167.4 ± 11.2 deg does not 
agree with either of the previously reported results 
lAnderson etail <f2010h found A = 210.7±g 1 9 5 deg, al- 
though this result was based o n only 3 data poin ts. 
More significant is the result by ITriaud et al.1 ()2010h of 
A = 212.7tgg deg, based on two data sets densely sam- 
pling the tra nsit. The diff e rence between our result for 
A and that of ITriaud et all (pOlol) is 45 ± 13 deg, or 3.5cr 
away from zero. 

We do not know the reason for the discrepancy. We 
have sparser covera ge of the spectroscopic orbit than did 
ITriaud et al.l 1)20101), but this did not bias our results: we 
used the results of lAnderson et all £2010) to set a prior 
for K, and also found that the results for K and A are un- 
correlated. We also checked to see if the gap in our time 
coverage when WASP- 17 was at zenith has biased our 
results: we fitted fake data with A = 210?7 and the same 
time coverage and velocity errors as our MIKE data, and 
found A = 206° ± 9°, i.e., the result was not biased. We 
also fitted the ITriaud et al.l (|2010l ) data ourselves and 
found A = 208° ± 12°. Although this reduces the discrep- 
ancy in results to 2.5er, it indicates that the differences 

5 For e ase of comparison, w e added 360 deg to the result 
quoted by Anderson et al. (2010), and we converted the result of 
ITriaud et all pOlOh onto our coordinate system using A = — j3. 



in A values cannot be attributed solely to differences in 
fitting procedures. An important step in our data anal- 
ysis was fitting and removing the airmass trend in our 
radial velocity data. We fitted for the airmass terms si- 
multaneously with the Rossiter-McLaughlin model, and 
therefore our results do take into account any correlations 
between the airmass parameters and the other parame- 
ters including A. However, these errors are internal to the 
choice of our model, which is linear in airmass and or- 
der number. This model seems reasonable and provides 
a satisfactory fit to the data; however it is impossible 
to exclude the possibility that the true dependence is 
more subtle and that this has altered the shape of the 
Rossiter-McLaughlin effect. 

According to the ephemeris of lAnderson et ail ([2010), 
the predicted time of conjunction for the event we ob- 
served is HJD 2,455,328.6814 ± 0.0017, which is ear- 
lier by 15 ± 5 min than the time we measured. Either 
the transits are nonperiodic, or the uncertainties in at 
least some of the transit times were underestimated. A 
st raight-line fit to the epo chs and transit times in Table 5 
of lAnderson et al.l (|2010D gives x = 36 with 11 degrees 
of freedom, i.e., statistically different than random. For 
this reason, it is likely that the true uncertainty in the 
period is larger than the previously reported uncertainty, 
by a factor of about ^/36/ll or 1.8. 

Consequently, we conclude that our t iming offset from 
the ephemeris of lAnderson et all ()2010D is at most a 1.5c 
discrepancy. The uncertainty in the time of conjunc- 
tion would be best addressed in the future by obtain- 
ing high-quality photometric data simultaneously with 
spectroscopic observations. However we note that this 
ephemeris discrepancy cannot acco unt for the differenc e 
between our result for A and that of lTriaud et al.l (2010). 

Our result for A is consistent with a retrograde orbit, 
but it must be remembered that A is a sky-projected 
quantity. The true angle tp between the vectors is given 

by 

cos ijj — cos i+ cos i + sin i± sin i cos A, (4) 

where i and i+ are the line-of-sight inclinations of the 
orbital and stellar angular momentum vectors, respec- 
tively. Using i — 86.95lQ;g3 deg from lAnderson et al.l 
(2010), and supposing to be drawn from an "isotropic" 
distribution (uniform in cosi*), we find ip > 92?6 with 
99.73% confidence. In this sense, the WASP-17b orbit is 
very likely to be retrograde (ip > 90°), although nearly 
pole-on configurations (tp ~ 90°) are possible. 

The assumption of an isotropic distribution in i± ne- 
glects our prior knowledge that main-sequence stars have 
somewhat predictable rotation rates, and therefore that 
the measured valu e of vsinij, also b ears information 
about i+ . Recently, Schlaufman (201Q|) used this insight 
to seek evidence for small values of sin (and therefore 
large spin-orbit misalignments along the line of sight) 
among all the transit hosts. He did not identify WASP- 
17 as one of the systems with a likely small value of sin i+. 
In light of this analysis, our observation of a large A fa- 
vors a more nearly retrograde orbits over polar orbits. 

4. DISCUSSION 

It is now ten years since the first reported detection of 
the Rossiter-McLau ghlin effect due to a transiting planet 
(|Queloz et al.l 12000). As with many other properties of 



4 



Bayliss et al. 



Fig. 3. — Results for the model parameters, based on our 
MCMC analysis of the radial velocity data, for parameters which 
show a degree of correlation with A. The contours represent 68%, 
95%, and 99.73% confidence limits. The quantity AT C is defined 
as T c minus the optimized value of HJD 2,455,328.6918. 

extrasolar planets, the Rossiter-McLaughlin effect obser- 
vations have thrown up surprises, with systems that ap- 
pear very different from those in our own Solar System. 
The Hot Jupiter HAT-P-7b was found to be on a po- 
la r or retrogra de orbit, b y independe nt measurements 
of lNarita et ail (|2009l ) and I Winn etaLl (|2009l ). Similarly 
WASP-17b has now been confirm ed to be in a retrog rade 
orbit by two independent groups (|Triaud et al.ll201Cl this 
work). It now seems that at least a fraction of Hot 
Jupiters are migrating by mechanisms other than disk 
migration, or that proto-planetary disks are frequently 

misaligned wit h their stars. 

The work of iFabrvckv fc Winn! (|2009h suggested that 
we may be seeing two different populations of Hot 
Jupiters that have migrated via different mechanisms. 
The results of this work adds weight to that suggestion, 
however a larger statistical sample will be needed before 
this can be robustly confirmed. Only a larger sample 
of planets with known A values will allow correlations 
with parameters such as planetary mass, planetary ra- 
dius, metallicity, and stellar mass to be robustly tested. 



W inn et all (|2010l ) have proposed that misaligned Hot 
Jupiters occur preferentially around hot stars (T c g > 
6250 K), of which WASP-17 is a supporting example. 
Fortunately the discovery rate of transiting planets or- 
biting bright stars is set to rise du e to large-scale g r ound - 
based surveys su ch as Super WASP Cameron et al.l (|2007l ) 
and HAT-South (|Bakos et al . 2009), and proposed s pace- 
based pr ojects such as TESS (iDeming et al.l l2009t ) and 
PLATO (|Catala fc PLATO Consortiumll2008lh 
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TABLE 1 

WASP-17 Radial Velocities (order averaged) 



Time 


Radial 


Radial velocity 


HJD 


velocity (m s _1 ) 


uncertainty (m s _1 ) 


2455328.57331 


-9.88 


38.75 


2455328.58273 


-4.73 


32.47 


2455328.59219 


9.70 


27.08 


2455328.60194 


26.88 


27.51 


2455328.61135 


-49.67 


37.43 


2455328.62035 


-87.48 


24.96 


2455328.62971 


-68.18 


20.08 


2455328.63921 


-110.63 


40.87 


2455328.64876 


-141.71 


30.59 


2455328.65817 


-124.57 


35.23 


2455328.66761 


-92.68 


27.30 


2455328.67700 


-135.62 


27.72 


2455328.68633 


-42.40 


34.15 


2455328.69570 


-13.34 


24.05 


2455328.70557 


50.30 


23.05 


2455328.71481 


37.16 


15.00 


2455328.74077 


164.59 


17.62 


2455328.75012 


86.38 


23.37 


2455328.75939 


86.77 


38.93 


2455328.76863 


31.27 


30.88 


2455328.77785 


6.61 


39.38 


2455328.78711 


—32.55 


18.04 


2455328.79641 


-30.80 


24.08 


2455328.80593 


-15.31 


24.74 


2455328.81546 


-23.01 


32.34 


2455328.82536 


-25.39 


32.21 


2455328.83495 


-18.66 


20.94 


2455328.84430 


6.76 


26.41 


2455328.85358 


50.02 


32.97 


2455328.86289 


-15.15 


25.22 


2455328.87214 


-45.57 


30.49 


2455328.88145 


13.95 


15.26 


2455328.89069 


-44.50 


24.20 
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TABLE 2 
Model Parameters for WASP- 17b 



Parameter Value 



Parameters determined mainly by the 
new radial velocity data 



Projected spin-orbit angle, A [dcgj 


167.6 ± 11.4 


Projected stellar rotation rate, vsini* [km s -1 ] 


8.61 ± 0.45 


Radial velocity offset, cq [m s -1 ] 


-0.267 ±0.074 


Order dependence of radial velocity offset, ci [m s _1 order -1 ] 


-0.79 ±0.11 


Airmass term, C2 [m s~ 1 airmass -1 ] 


0.233 ± 0.061 


Order dependence of airmass term, C3 [m s — 1 airmass - 1 order -1 ] 


0.738 ± 0.093 


Parameters determined mainly by priors 


Velocity semi-amplitude, K [m s 1 ] 


55.9 ± 5.0 


Time of conjunction, T c [HJD-2,450,000] 


5, 328.6919 ±0.0027 


Transit depth, (R p /R+) 2 


0.01676 ± 0.00027 


Total transit duration, tiv — t\ [hr] 


4.368 ± 0.037 


Partial transit duration, tiv — fin [hr] 


0.566 ± 0.039 



